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SUMMARY 

A detailed benzene oxidation mechanism, which uses the qualitative paths 
outlined by previous Investigators, Is presented. Ignition delay times for 
argon-diluted mixtures computed with this mechanism are In satisfactory agree- 
ment with experimental values obtained from pressure-time profiles for a wide 
range of Initial conditions. An experimental temperature versus time profile 
measured In a nitrogen diluted mixture has also been successfully computed. 
Additional computations have qualitatively matched several experimental species 
concentration versus time profiles for the nitrogen diluted reaction. Applica- 
tion of sensitivity analysis allowed the approximate determination of the rate 
constants for two Important reactions: 


C 6 H 5 + °2 


W ♦ 0 


k 7 . 4.5X10 12 e-' 5000 cal/RT c» 3 mor's " 1 


and 


c 5 H 5 ♦ 0 2 



c 5 h 5 o ♦ 0 


k 


10 


2.0xl0 12 


-20000 cal/RT 3 .-l c -l 
e cm mol S 


INTRODUCTION 


The Increasing Importance of aromatics in the practical hydrocarbon fuels 
of today has accelerated research aimed at Improving our understanding of the 
high-temperature oxidation of these compounds. This knowledge Is important in 
controlling the combustion and emission characteristics of both advanced air- 
craft propulsion and future ground-based gas turbines. The oxidation and 
pyrolysis of the simplest aromatic, benzene, have been studied by several 
Investigators (refs. 1 to 4). A recent review paper (ref. 5) on aromatic hydro- 
carbon oxidation presents experimentally measured temperature and composition 
profiles for oxidation of one benzene-oxygen-nitrogen mixture in a highly tur- 
bulent flow reactor. In some recent work (ref. 6) Ignition delay time measure- 
ments were reported for benzene-oxygen-argon mixtures Ignited behind a 
reflected shock. The experimental conditions covered a wide range of Initial 
composition, temperature, and pressure. The ignition delay times were deter- 
mined from pressure versus time profiles. 



The purpose of the present paper is to develop a detailed mechanism which 
reproduces the results of references 5 and 6 over the widest possible range of 
initial conditions. Computations were performed with the new NASA Lewis 
general chemical kinetics and sensitivity analysis code (refs. 7 to 9). To 
achieve the best possible agreement sensitivity analysis was used to identify 
the most important reactions. Adjustments were then made to some of the 
reaction rate constants in the mechanism. However, only rate constants with 
large uncertainties were changed significantly. These adjustments were made 
mainly to reactions involving benzene and its fragments. The rate constants 
for reactions involving other hydrocarbons, e.g., acetylene and methane, were 
used at their literature values. In the sections which follow, we present 
comparisons of computed and experimental results and discuss the mechanism 
finally obtained for optimal matching with the latter results. We also 
describe the sensitivity study which determined the reactions which have the 
greatest effect on the computed results. 


BENZENE OXIDATION MECHANISM 

We have followed the general scheme outlined qualitatively in references 1 

and 5. The oxidation process can be divided into three stages: 

Initiation: The formation of oxygen and hydrogen atoms and OH radicals as 
well as some phenyl radicals by the pyrolysis of benzene and its reaction 
with molecular oxygen. 

Induction: Rapid formation of phenyl radical and Its oxidation and pyrolysis 

to form primarily phenoxy radical, CgHgO; further reaction of C 5 H 5 O 
to form phenol, cyclopentadi ene , and acetylene. 

Combustion: Oxidation of the stable intermediate hydrocarbons and radical 

recombinations which cause the main heat release. 

In the initiation region the following reactions were assumed: 


i * °2 ** C 6 H 5 * H0 2 

( 1 ) 

C 6 H 6 % C 6 H 5 ♦ H 

( 2 ) 

k 3 


H + 0 2 OH + 0 

(3) 


In the induction region the reactions used were as follows: 

H + C 6 H 6 t=^ C 6 H 5 + H 2 
0 + C g Hg ^ CgHg + OH 


(4) 

(5) 


2 



OH + C 6 H 6 


c 6 h 5 ♦ h 2 o 


( 6 ) 


We have written reactions 4 through 6 as simple abstraction processes, although 
the possibility of addition processes has been suggested (ref. 5). 


C 6 H 5 +0 2 ^ C 6 H 5 0+0 


C 6 H 5 0 + M Ji C 5 H 5 + CO + m 


cyclo- 

pentadienyl 


CO + OH C0 2 + H 


C 5 H 5 + °2 C 5 H 5° + 0 


C 5 H 5 0 + M JJ C 4 H 5 + CO + M 
butadienyl 


(7) 

( 8 ) 


(9) 

(10) 

( 11 ) 


c 6 h 5 oh 



c 6 h 5 o + H 


(12) 


C 5 H 5 + C 6 H 6 ^ C 5 H 6 + C 6 H 5 (13) 


cyclo- 


pentadiene 


^14 

oh + c 6 h 5 oh c 6 h 5 o + h 2 o 

(14) 

c 6 h 5 c 4 h 3 , c 2 H 2 

(15) 

k lfi 

H + C 4 H 3 C 4 H 2 + H + M 

(16) 

^17 

C 5 H 6 + 0 + C 5 H 5 0 + H 

(17) 

^1 8 

c 4« 5 C 2 H 3 , C 2 H 2 

(18) 


These are the significant steps which occur in the initiation and induction 
regions. Additional reactions, discussed later, are needed to model the com- 
bustion region and determine the concentrations of intermediate species. The 
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complete set of reactions Is necessary for modeling the temperature and com- 
position profiles. 


REACTION RATE CONSTANTS 

In this section we discuss the selection of rate constants for those 
reactions Involving benzene and Its decomposition products which have been 
studied In recent experimental work. Sources for rate constant of combustion- 
zone reactions are also given. 


Benzene Pyrolysis Reaction 

The pyrolysis of benzene has been studied recently In two Investigations, 
by Hsu, Lin, and Lin (Ref. 4) and by Kiefer et al. (Ref. 10). This work has 
shown that reaction (2) Is the only Important path at high temperature. The 
pyrolysis rate constant was measured as a function of temperature and pressure. 
For the temperature and pressure range of the present computations we have used 
the expression given In reference 4, which agrees well with results of 
reference 10 . 


Radical Attacks on Benzene 

Reactions 4, 5 and 6 are the main reactions for generating phenyl radical. 
We have used the rate constant of reference 10 for the H + reaction. 

The recently reported experimental values of references 11 and 12 were used for 
the 0 + C 5 H 5 and OH + C 5 H 5 reactions. In a later section we discuss sensitiv- 
ity analysis computations which showed the pressure and temperature profiles to 
be quite Insensitive to moderate variations of these three rate constants from 
their literature values. This Is true, even though these reactions are all 
Important sources of phenyl radical. 


Other Cfc and C 5 Reactions 

The rate constants of some reactions of phenyl and phenoxy (C 5 H 5 O) radical 
have been measured recently. The previously quoted work of reference 10 has 
determined rate constant expressions for the phenyl radical pyrolysis, reaction 

15. Their high-pressure unlmolecular reaction expression was used In the pre- 
sent computations. The estimated rate constant for pyrolysis of C 4 H 3 , reaction 

16, Is also taken from reference 10. The rate constant for the phenoxy pyrol- 
ysis, reaction 8 , Is the value of Lin and Lin (ref. 13). The rate constants 
for other reactions In the Induction region had to be either estimated or 
adjusted to match the experimental data. The sensitivity analysis computations 
for accomplishing this task are described below. 
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reference 14. The reactions of acetylene, which are important in the combus- 
tion region, were all taken from reference 15. The complete mechanism con- 
sisted of 110 reactions among 35 species. All reactions and rate constants 
used for these computations are given in table I. 


SENSITIVITY ANALYSIS STUDY 

To identify the important reactions which control the oxidation of benzene 
and its fragments in the induction zone, an extensive sensitivity analysis was 
performed. Normalized sensitivity coefficients were computed using the 
decoupled direct method described in references 8 and 9. These coefficients 
give the variation of all species concentrations, temperature and pressure with 
changes in the individual rate constant parameters Aj, n j , and Ej of the 

n 1 

modified Arrhenius equation k^ = A^T J exp (-E^/RT) for reaction j. In the 

present work sensitivity coefficients of several species concentrations plus 
temperature and pressure with respect to the Aj parameters were used to 
determine the important reactions. Table II shows these coefficients for the 
nitrogen diluted stoichiometric oxidation of benzene (initial temperature = 

1115 K) in the turbulent flow reactor described in reference 5. The reaction 
is at constant pressure so no pressure sensitivity coefficient is given. The 
results show that the reactions of C 5 H 5 and C 5 H 5 radicals with molecular 
oxygen have greatest effect on temperature and several species composition 
profiles. Table III shows sensitivity coefficients for a typical stoichio- 
metric benzene oxygen-argon shock ignition (initial temperature = 1405 <) from 
reference 6 . Sensitivity coefficients of pressure are now included. Again the 
dominant reactions are the CgHs and C 5 H 5 oxidations. Computations were also 
performed for a second shock ignition at equivalence ratio = 2 and an initial 
temperature of 1600 K. Sensitivity coefficients for four dependent variables 
and five reactions are shown for both conditions in table IV for ease of com- 
parison. The data in table IV show that computed results should be most 
sensitive to the rate constant of reaction 7, the phenyl oxidation, at both 
conditions. However, there is a significant change in the sensitivity coef- 
ficients for reaction 15, the phenyl pyrolysis. These coefficients change 
from very small values for the stoichiometric lower temperature reaction to 
much larger values for the rich higher temperature reaction. The effect of 
reaction 15 on the pressure profile is now Indicated to be greater than that 
of reaction 10 . 

The predictions of the sensitivity coefficients were tested by using the 
"brute force" or indirect method of sensitivity analysis. The preexponential 
factor for each reaction in table IV was changed by factors of 2 and 0.5 and 
the ignition computation was repeated for each change. The effect of these 
rate constant variations on the computed pressure profile is shown in table V. 
Given in this table are values of xp, the ignition-delay time measured from 
the pressure versus time profile as described in detail in the next section. 

For each reaction we show the xp values for k^p, the rate constant listed 
in table I, and also xp for 2k s td and k s t<j/2. Also shown are the percent 
changes in x p for each rate constant variation, the computed results in 
table V confirm all the sensitivity predictions of table IV. In particular, 
changing k -| 5 for the lower-temperature stoichiometric mixture has a very 
small effect on x p . The same change in k ]5 has a very significant effect 
on xp for the high-temperature rich mixture. 
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The results of these sensitivity computations were used in the following 
way to adjust the rate constant parameters. The Aj and Ej values for 
reactions 7 and 10 were adjusted in several iterations to give the best match 
to the ignition delay times of reference 6 and the temperature versus time 
profile of reference 5 for the benzene-oxygen-ni trogen reaction. It was found, 
from the matching of ignition delay times, that the activation energy of the 
dominant reaction 7 should be as high as possible. This reaction is an exo- 
thermic process and we have used the upper limit, suggested in reference 1 , of 
15 kcal/mole. The C 5 H 5 oxidation, reaction 10, is about 10 kcal/mole 
endothermic at 1200 K. An activation energy of 20 kcal/mole was found to give 
the best shape to the computed temperature versus time profile for the benzene- 
oxygen-nitrogen reaction. The prexponential factors of both these reactions 
were then adjusted to give the best match to both the ignition-delay times and 
the temperature versus time profile. Table V and the sensitivity coefficients 
of table II show that reaction 7 strongly controls the t p computation and 
that reactions 7 and 10 both have Important effects on the temperature versus 
time profile. Therefore the prexponential factor A 7 was determined by match- 
ing the ignition-delay times and A] g was then adjusted to give the best match 
to the temperature versus time profile. The values obtained are: 

k 7 = 4.5xl0 12 e - 15000 cal/RT cm 3 moHS ' 1 

k^o = 2.0xl0 12 e ~ 20000 cal/RT cm 3 mol - 1 S -1 

Other rate constants for reactions involving cyclopentadiene, phenol and 
their fragments were estimated and adjusted to give the best agreement with 
reported experimental phenol and cyclopentadiene profiles in reference 5. 


DESCRIPTION OF COMPUTATIONAL PROCEDURE 

In this section we describe briefly the apparatus and experimental 
procedures reported in references 5 and 6 . We then describe our mathematical 
model used to simulate each of these reacting systems. 


Turbulent Flow Reactor 

Temperature and composition versus time profiles for the reaction of a 
stoichiometric benzene-oxygen mixture in a nitrogen atmosphere were obtained 
in the Princeton University flow reactor (ref. 5). A detailed description of 
this turbulent, chemical-rate limited device is given in reference 16. The 
reactants were highly diluted In nitrogen to eliminate diffusion effects. 
Profiles of temperature and composition with respect to distance were converted 
to time profiles by taking into account the flow velocities within the reactor. 
The exact zero of reaction time is arbitrary and was taken as point of fuel 
injection into the hot oxidant stream. Reference 16 indicates that the reactor 
is operated at one atmosphere pressure. In the computations, therefore, the 
flow reactor was modeled as a constant pressure homogeneous batch reaction. 

For all kinetic computations of this work, the thermodynamic data are from the 
NASA Lewis Research Center data base which is part of the NASA Chemical Equili- 
brium Composition Computer code of Gordon and McBride (ref. 22). The recently 
reported data of Burcat, Zeleznik and McBride (ref. 23) were used for phenyl 
and phenoxy radicals. New thermodynami c data for several benzene pyrolysis and 

6 



oxidation fragments (C4H3, C4H5, C5H5 and C5H5O) were estimated by 
Bonnie J. McBride. Computed time profiles of temperature and several species 
concentrations were compared to the profiles presented In reference 5. 


Shock-Tube Ignition Experiments 

The Ignition delay times of Ref. 6 were obtained from measured pressure 
versus time traces following Ignition behind a reflected shock wave. The time 
Interval between Ignition and the first observed "significant" pressure rise 
was taken as the Ignition delay time, x. To match the experimental technique 
as closely as possible the computed Ignition delay time, x p , was obtained from 
the computed pressure versus time curve. The kinetic computations were per- 
formed assuming a constant-volume batch reaction behind the reflected shock. 
Before performing the kinetic computations the Initial reflected shock condi- 
tions In reference 6 were recomputed. A small correction for attenuation of 
the reflected shock velocity was applied to each data point. The reflected 
shock Initial tempertures and pressures (T5 and P5) were then recomputed 
using the Chemical Equilibrium Code of reference 22. Complete time profiles of 
temperature, pressure, and all species concentrations were then computed. A 
typical computed pressure versus' time curve Is shown In figure 1. Also shown 
Is the method used to determine the x p value by the Intersection of two 
extrapolated lines. This method approximated as closely as possible the exper- 
imental technique. Only 35 of the experimental points reported In reference 6 
were used for our comparisons. The other points were excluded from consider- 
ation for one of two different reasons. First, all points for which x was 
less than 100 psec were excluded. These values are Inaccurate because of pres- 
sure disturbances which propagate In the shock tube ahd. cause nonuniform heat- 
ing of the gas behind the reflected shock. Details of this effect are given In 
reference 19. Second, all the data taken for the experimental mixture with 
equivalence ratio of 0.25 were excluded because It was very difficult to deter- 
mine accurate x values from these pressure traces. The pressure rise was 
very poorly defined for all these very weak Ignitions. 

Table VI lists the five mixtures from reference 6 which were used In the 
present comparisons along with the Initial temperature range of the reflected 
shocks. Mixtures 3 and 4 are Identical for all practical purposes and are 
treated as one condition, equivalence ratio 1.0 with about 85 percent argon 
dilution. Note that comparison of the results for mixture 2 with those for 
mixtures 3 and 4 gives the effect of simply diluting the mixture with argon at 
a constant equivalence ratio. Initial pressures for the reflected shocks 
ranged from 1.7 to 7.1 atm. 


RESULTS AND DISCUSSION 

Comparison of Computed and Experimental Ignition Delay Times 

Plots of computed and experimental Ignition delay time versus reciprocal 
of temperature are shown In figures 2 to 6 for the four different initial con- 
ditions given in table VI. The Individual computed and experimental points 
are shown as well as least-squares lines for each set of points. They are 
fitted to the empirical equation 
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X 


Ae* E/RT 


(1) 


or 


log T a 1og A + 2.30 T R T (2) 

where R Is the universal gas constant and AE Is an activation energy term 
which measures the temperature dependence of t for a fixed set of initial 
concentrations. A detailed comparison of experimental and computed results is 
shown in table VII, which also Includes an error analysis. The percent differ- 
ence between each computed and experimental value Is given along with the 
percent standard deviation for each set of initial concentrations used. This 
standard deviation Is defined by 


x 100 % <3) 

where N Is the number of points. A value of a for all 35 points used Is 
also given. 

Examination of this table shows satisfactory agreement between the com- 
puted and experimental ignition delays. This Is especially true when one con- 
siders the uncertainty In some of the experimental measurements. For example, 
table VII shows three Instances In which the experimental results do not show 
the expected decrease of ignition delay time with Increasing temperature. The 
agreement between experiment and computation is poorest for the lean mixture 
(<J> * 0.5) and steadily Improves as $ increases to 2.0. Similar behavior 
Is also observed for the temperature dependence of Ignition delay. This can 
be seen from table VIII which lists the values of AE computed from equation 
(1) for both the experimental and computed results. The percent deviation for 
the activation energies decreases steadily as equivalence ratio changes from 
0.5 to 2.0. From these comparisons It appears that one or more reactions are 
missing from the mechanism. These reactions are much more important in the 
lean mixtures than In the rich ones. However, It should also be pointed out 
that the experimental Ignition delay times for the lean mixtures are more 
uncertain than those for the rich mixtures. From figures 2 and 5 It can be 
seen that the computations tend to overpredict the ignition delay time at 
4> - 0.5 and underpredict It at 4> => 2.0. 

Figure 6 shows computed and experimental Ignition delay plots for two 
mixtures with equivalence ratio of 1.0 but with different amounts of argon 
diluent. The computed results show a much smaller effect of inert dilution 
than observed experimentally. The effect of Inert gas dilution could not be 
changed by variation of either the A factor or the activation energy of any 
uncertain rate constants. 
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Comparison of Computed and Experimental Temperature and Composition Profiles 

Figures 7 to 12 show computed temperature versus time and several composi- 
tion versus time profiles along with experimental data for the nitrogen-diluted 
stoichiometric CgHg-C^ reaction reported in reference 5. The initial mixture 
temperature was 1115 K and the pressure was 1 atmosphere. The computed and 
experimental temperature profiles in figure 7 are in excellent quantitative 
agreement. The computed cyclopentadiene profile of figure 8 is in semiquan- 
titative agreement with the experimental profile. Both the peak concentration 
of C 5 H 5 and the time at which it is reached are predicted quite accurately. 

The computed profiles for phenol (fig. 9) benzene (fig. 10) and carbon monoxide 
(fig. 11 ) agree only with the qualitative trends of the corresponding exper- 
imental profiles. For phenol, the computation matches the magnitude on the 
peak concentration, but this value occurs at a later reaction time. The mecha- 
nism predicts a very sudden disappearance of benzene, instead of the gradual 
decay observed experimentally. The experimental rise of carbon monoxide con- 
centration and its peak value are greater than the corresponding computed 
quantities. There is semi-quantiati ve agreement between the computed and 
experimental carbon dioxide profiles (fig. 12). The computed profile is 
approximately linear whereas the experimental profile shows a long induction 
period. However, the same steady-state value found experimentally is given by 
the computed profiles. 

An acetylene versus time profile is also given in reference 5 for this 
benzene oxidation experiment. The predicted acetylene peak concentration is 
about one order of magnitude lower than the reported experimental value. The 
sensitivity coefficients of table II show that this concentration is controlled 
by reactions 7 and 10, as are many other concentrations . No reasonable changes 
in any of the reactions involving acetylene changed its profile significantly. 
This behavior is different from that for cyclopentadiene and phenol. Although 
these latter two concentrations are also strongly controlled by reactions 7 and 
10 , they were significantly improved by modest adjustments of the rate con- 
stants for reactions 12 and 13, which are quite uncertain. 


CONCLUDING REMARKS 

This work has presented a detailed quantitative mechanism for the oxida- 
tion of benzene. The mechanism satisfactorily computes experimental results 
for both argon and nitrogen diluted systems. Fair to good agreement was 
obtained between computed and experimental ignition delay times measured over 
a wide range of temperature, pressure and composition for shock-heated benzene- 
oxygen-argon mixtures. The mechanism also does a fair job of computing meas- 
ured temperature and composition profiles for a nitrogen-diluted benzene 
oxidation in a quite different system, namely a turbulent flow reactor. 

The present work has verified the reaction paths outlined in references 1 
and 5 and presents rate constant expressions for two important reactions. 
Sensitivity analysis has identified these reactions, the oxidations of C 5 H 5 
and C 5 H 5 radicals, as being dominant in the ignition of benzene-oxygen mix- 
tures. The expressions for these reaction rate constants have to be considered 
approximate because the mechanism is still incomplete. However, the mechanism 
predicted one temperature versus time profile accurately, in addition to the 
satisfactory ignition delay time predictions. Therefore, it could be used for 
heat-release rate predictions in practical, well-mixed, benzene-air combustion 
systems. 
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TABLE I. - BENZENE OXIDATION MECHANISM 


No. 

Reaction 

BB 

n 

E a 

calories 
per mole 

Source 

1 

C 6 H 6 + °2 ■ C 6 H 5 * H0 2 

6.31X10 13 

0 

60 

000 

Ref. 

2 

B 

C 6 H 6 ■ C 6 H 5 + H 

5-OxlO 15 

0 

108 

000 

Ref. 

4 

B 

H + 0 2 = OH + 0 

1 .26xl0 14 

0 

16 

300 

Ref. 

14 

■ 

C 6 H 6 * H ' C 6 H 5 * H 2 

1 4 

2.5x10 

0 

16 

000 

Ref. 

10 

1 

C 6 H 6 * 0 = C 6 H 5 ♦ OH 

2 .78x1 0^ 3 

0 

4 

910 

Ref. 

11 

6 

C 6 H 6 ♦ OH . C 6 H 5 ♦ H 2 0 

1 3 

2.13x10 

0 

4 

580 

Ref. 

12 

7 

C 6 H 5 * °2 - C 6 H 5° f 0 

12 

4.5x10 

0 

15 

000 

This 

work 

8 

c 6 h 5 o - c 5 h 5 ► CO 

2.5X10 11 

0 

43 

900 

Ref. 

13 

9 

CO + OH = C0 2 ^ H 

4.17X10 11 

0 

1 

000 

Ref. 

14 

m 

C 5 H 5 *■ 0 2 = C 5 H 5 0 * 0 

1 2 

2.0x10 

0 

20 

000 

This 

work 

11 

C 5 H 5 0 ♦ M - C 4 H 5 + CO t M 

7.59X10 13 

0 

15 

000 

Est. 


12 

C 6 H 5 0H = C 6 H 5 0 * H 

1 3 

6.0x10 

0 

1 88 

000 

Est. 


13 

C 5 H 5 + C 6 H 6 = C 5 H 6 + C 6 H 5 

2 .oxTo” 

1 0 

10 

000 

Est. 


14 

oh t c 6 h 5 oh = c 6 h 5 o + h 2 o 

1 2 

8.0x10 

0 

5 

000 

Est. 


15 

C 6 H 5 * C 4 H 3 + C 2 H 2 

T ,58xl0 15 

0 

82 

000 

! Ref. 

10 

16 

M f C^Hj - C 4 H 2 + H ♦ M 

3.3xl0 51 


63 

000 

Ref. 

10 

17 

C 5 H 6 «■ 0 = CgH^O * H 

1 2 

5.0x10 

0 

10 

000 

Est. 


18 

c 4 h 5 - c 2 h 3 ♦ c 2 h 2 

1 3 

1 .4x10 

0 

32 

900 

Est. 


19 

C 4 H 2 *• 0 * C 2 HO + C ? H 

1 3 

1.0x10 

0 


0 

Ref. 

2 

m 

C 4 H 2 f 0 - CO + C 3 H 2 

1 2 

1.2x10 

0 


0 

Ref. 

15 


12 













TABLE I. - Continued. 


No. 

Reaction 

A 

cm 3 , mole, sec 

n 

Ea 

calories 

Source 

21 

C 4 H 2 + OH + HCO + CjH 2 

3.0x10 13 



Ref. 

15 

22 

C 2 H 4 + M + C 2 H 2 ♦ H 2 ♦ M 

9 . 33xT 0 1 6 



Ref. 

15 

23 

C 2 H 4 + OH = C 2 H 3 + H 2 0 

1 2 

4.79x10 

0 

1 230 

Ref. 

18 

24 

C 2 H 4 + OH + CH 3 + CH 2 0 

2. 00x1 0 1 2 

0 

960 



25 

C,H. + 0 + CH„ + HCO 
c 4 3 

1 2 

3.31x10 

0 




26 

C 2 H 4 + 0 = CH 2 0 ♦ CH 2 

2.51xl0 13 

0 




27 

C 2 H 3 + H = C 2 H 2 + H + H 

14 

7.94x10 

0 




28 

C 2 H 3 * 0 2 = C 2 H 2 + H0 2 

1 . 58x1 0 1 3 

0 


Ref. 

15 

29 

C 2 Hj + H . C 2 H 2 + H 2 

1 2 

6.0x10 c 

0 




30 

C 2 H 3 + OH = C 2 H 2 + H 2 0 

5.00xl0 12 

0 




31 

c 2 h 3 + CH 2 = c 2 h 2 ♦ CH 3 

3.00X10 13 

0 




32 

C 2 Hj + C 2 H = C 2 H 2 + C 2 H 2 

3.00X10 13 

0 

0 



33 

C 2 Hj + 0 = C 2 H 2 0 + H 

3.30X10 13 

0 

0 



34 

CH 2 + CH 2 = C 2 H 2 ♦ H 2 

4.00X10 13 

0 

0 



35 

CH 2 + CH 2 = C 2 H 3 + H 

5.01X10 12 

0 


Ref. 

18 

36 

CH 2 t OH = CH + H 2 0 

2.51X10 11 

.67 


Ref. 

15 

37 

CH 2 t 0 = CH + OH 

2-OOxlO 11 

.68 




38 

CH 2 + 0 2 = C0 2 + 2H 

1.59X10 12 

0 




39 

C 2 H 2 + M = C 2 H + H + M 

4.17X10 16 

0 




40 

C 2 H 2 + C 2 H 2 = C 4 H 3 + H 

2.00x1 0 1 2 

0 




41 

C 2 H 2 + 0 2 = C 2 H0 + OH 

2.00x10® 

1.5 




42 

C 2 H 2 + 0 = C 2 H t OH 

3.1 6x1 0 1 ® 

-.6 




43 

c 2 h 2 + o = ch 2 ♦ CO 

2.20x10*® 

1.0 

2 580 



44 

C 2 H 2 + 0 = C 2 H0 + H 

3.55x1 0 4 

2.7 




45 

C 2 H 2 + OH = C 2 H + H 2 0 

6.31xl0 12 

0 


Ref. 18 
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TABLE 



Reaction 


C 2 H 2 + OH . C 2 H 2 0 4 H 


C 2 H 2 + C ? H . C 4 H 2 * H 


C 2 H 2 + CH ? = C 3 H 3 + H 


CjH 3 4 H 4 M . CjH 4 4 M 

50 

c 2 h 2 o * oh = ch 2 o + HCO 

51 

C 2 H 2 0 + OH = C 2 HO + H 2 0 

52 

C 2 H 2 0 + h = ch 3 + CO 

53 

C 2 H 2 0 4 H . C 2 H0 + H 2 

54 

C 2 H 2 0 + 0 = C 2 H0 + OH 

55 

C 2 H 2 0 + 0 = CH 2 0 + CO 

56 

C 2 H 2 0 + m = ch 2 <• CO + m 

57 

C-jHO + 0 2 = 2C0 + OH 

58 

C 2 HO + 0 » 2C0 4 H 

59 

CjHO + OH . 2HC0 

60 

c 2 ho <■ h = ch 2 + CO 

61 

c 2 ho + ch 2 . c 2 h 3 * CO 

62 

C 2 H0 + CH 2 = CH 2 0 4 C 2 H 

63 

C 2 HO 4 C 2 HO = C 2 H 2 + 2C0 

64 

C 2 H * OH . C 2 H0 + H 

65 

C 2 H + 0 2 = C 2 HO 4- 0 

66 

C 2 H + 0 = CO 4 CH 

67 

C 2 H 4 H 2 = C 2 H 2 4 H 

68 

CH 4 4 M = CH 3 4 H 4 M 

69 

ch 4 4 o 2 = ch 3 4 ho 2 

70 

CH 4 4 H = CH 3 4 H 2 


I 


Continued . 


A 

cm 3 , mole, sec 

n 

Ea 

calories 

Source 

3. 20x1 O 11 

0 

200 

Ref. 

15 

3. OOxlO 13 

0 

0 



1 .OOxlO 12 

0 

0 



2. OOxlO 13 

0 

0 



2.80xl0 13 

0 

0 



7 . 5x1 o' 2 

0 

3 000 



1 . 1 3x1 0 1 3 

0 

3 428 

Ref. 

18 

7.50X10 13 

0 

8 000 

Ref. 

15 

5. OOxlO 13 

0 

8 000 



2. OOxlO 13 

0 

0 



2. OOxlO 16 

0 

60 000 



1 2 

1.46x10 

0 

2 500 

\ 


12 

1.20x10 

0 

0 

Ref. 

18 

1. OOxlO 1 3 

0 

0 

Ref. 

15 

5. OOxlO 13 

0 

0 



3. OOxlO 1 3 

0 

0 



1 .OOxlO 13 

0 

2 000 



1 .OOxlO 13 

0 

0 



2. OOxlO 1 3 

0 

0 



5. OOxlO 13 

0 

1 500 



5. OOxlO 13 

0 

0 



4 . 09x1 0 6 

2.39 

864 



2. OOxlO 17 

0 

88 000 

Ref. 18 

7 .94x1 O 1 3 

0 

56 000 

Ref. 18 

1 . 26x1 0 1 4 

0 

11 900 

Ref. 20 




TABLE I. - Continued. 


No. 

Reaction 

A 

cm 3 , mole, sec 

n 

Ea 

calories 

Source 

71 

CH 4 + OH = CH 3 + H 2 0 

2.5xl0 13 

0 

5 010 

Ref. 

19 

72 

CH 4 + 0 = CH 3 + OH 

1.9xl0 14 

0 

11 720 

Ref. 

19 

73 

CH 3 + OH = CH 2 0 + H 2 

3.98X10 12 

0 

0 

Ref. 

17 

74 

CH 3 + OH * CH 3 0 + H 

2.00xl0 16 

0 

27 410 



75 

CH 3 + 0 2 = CHjO + 0 

4 . 79x1 0 1 3 

0 

29 000 



76 

CH 3 + 0 = CH 2 0 + H 

14 

1 .29x10 

0 

2 000 



77 

CH 3 + CHj = C 2 H 4 + H 2 

1 .OOxlO 16 

0 

32 000 



78 

CH 3 0 + M = CH 2 0 + H + M 

5.01X10 13 

0 

21 000 



79 

CH 3 0 + 0 2 = CH 2 0 + H0 2 

1 .OOxlO 12 

0 

6 000 



80 

CH 3 0 + H = CH 2 0 + H 2 

2 . 00x1 0 1 3 

0 

0 



81 

CH 3 + CH 2 0 = CH 4 + HCO 

1 .OOxlO 10 

.5 

6 000 



82 

CH 3 + HCO = CH 4 + CO 

3.02x1 0 1 1 

.5 

0 



83 

CH 3 + H0 2 = CH 3 0 + OH 

2. OOxlO 13 

0 

0 



84 

CH 2 0 + M = HCO + H + M 

3.31X10 16 

0 

81 000 



85 

CH 2 0 + OH = HCO + H 2 0 

1 2 

7. 59x10 c 

0 

170 



86 

CH 2 0 + H = HCO + H 2 

3 . 31 xl O 1 4 

0 

10 500 



87 

CH 2 0 + 0 = HCO + OH 

5.01X10 13 

0 

4 600 


r 

88 

HCO + H0 2 = CH 2 0 + 0 2 

1 .OOxlO 14 

0 

3 000 

Ref. 

17 

89 

HCO +M=H+C0+M 

14 

2.94x10 

0 

1 5 570 

Ref. 

20 a 

90 

HCO + 0 2 = CO + H0 2 

3. 31 xlO 1 2 

0 

7 000 

Ref. 

17 

91 

HCO + OH = CO + H 2 0 

1 .OOxlO 14 

0 

0 



92 

HCO + H = CO + H 2 

2. OOxlO 14 

0 

0 



93 

HCO + 0 = CO + OH 

1. OOxlO 1 4 

0 

0 



94 

CH + 0 2 = HCO + 0 

1 .OOxlO 13 

0 

0 

Ref. 

17 


Computed from k for the recombination reaction and the equilibrium 
constant expression given In Ref. 18. 


15 



No. 

Reaction 

A 

n 

Ea 

Source 

; 


cm 3 , mole, sec 

i 

calories 

i 



95 CO + 0 + M = C0 2 + M 

96 CO + 0 2 = C0 2 + 0 

97 CO + H0 2 = C0 2 + OH 

98 Ot H ? 0 = OH + OH 

99 0 + H 2 = OH ^ H 

100 H 2 + 0 2 = OH +■ OH 

101 0 + H0 2 = OH + 0 2 

102 OH + H0 2 - H 2 0 + 0 2 

103 H + H0 2 = H ? 0 + 0 

104 H + H0 2 = OH + OH 

105 H 2 + OH == H 2 0 + H 

106 H + 0 2 f H = H0 2 + M 

107 H + OH + M = H 2 0 + M 

108 H+0+M=0H+M 


109 H+H+M=H 2 +M 
1 1 o 0+0+M=0„+M 


5.89x10 
2.51x10 
5.75x10 
6.76x10 
2.95x10 
2.51x10 
5.01x10 
5.01x10 
5.01x10 
2.51x10 
2.19x10 
1 .51x10 
1 .41x10 
1.00x10 
3.02x10 
1 .91x10 


Ref. 17 
47 690 Ref. 17 

22 930 Ref. 17 

18 360 Ref. 17 

Ref. 14 
Ref. 17 




Reaction 

number 

j 


Reaction 


Sensitivity coefficient of dependent variable n: 





C £ H c 

+ C c H c + HO, 

6 6 

2 6 5 ; 

C c H c 

^ C c H c 4- H 

6 6 

6 5 

H + ( 

] OH + 0 

C 6 H 6 

+ h^c 6 h 6 ♦ H 2 

C £ H £ 

+ 0 ^ C H c + OH 

6 6 

6 5 

C r H, 

+ OH ^ C r H c + HJ 

6 6 

6 5 2 

C £ H c 

+ 0„ ^ C c H c 0 + 0 

6 5 

2 6 5 

C £ H J 

7 ^ C c H c + CO 

6 5 

5 5 

CO + 

OH ^ C0 2 +■ H 

C H 

+ 0, ^ C c H c 0 + 0 

5 5 

2 5 5 

C,H_0H ^ C.H_0 + H 

6 5 

b 5 

C r H, 

+ C £ H c ^ C c H e + 1 

5 5 

6 6 5 6 

OH + 

C,H OH ^ C £ H r 0 + 


6 6 6 5 

C,H r 

^ C . H_ + C„H„ 

6 5 

4 3 2 2 

C c H c 

+ 0 C c H c 0 + H 

5 6 

5 5 

C-H„ 

+ o ^ ch 2 * CO 

2 2 



c 2 h 2 + 0 * 

H ♦ 0 2 + M' 


^HO + H 


c 6 h 6 
-0.178/ 
.8959 
- . 0530 

.0058 

-.8467 

-.8013 

-5.680 

.1169 

.3546 

.4904 

-.1854 

-.1399 

.0065 


C 2 H 2 


0.2050 -0.1169 0 . 0244 

.0610 .5429 .3042 

.2073 -.2415 .1247 


.0056 

-.3351 

- .4225 
-6.331 

.0804 

.1428 

.2212 

- .0996 
-.1176 
-.0045 

4x1 0~ 5 
.0083 
.0456 


.0030 
-.1150 
-.1949 
-3.333 
- .0435 
-.2367 
.2088 
.3744 
-.0772 
.1317 

-.0001 

.0157 

-.3350 


c 5 h 6 c 6 h 6 oh 


0.1129 

-.4134 

-.0685 

-.0027 

.1204 

.1847 

3.089 

-.1051 

.2001 

.4070 

.6343 

.0439 

-.4036 


Aj/n (dn/aAj) at time = 35 msec 


CO 


0.0359 
.0469 
- .0025 




16 



























17 


0607 -.0009 -.0010 



TABLE IV. - SENSITIVITY COEFFICIENTS FOR SEVERAL REACTIONS FOR TWO 
SHOCK IGNITIONS 


Reaction 

number 

Reaction 

Sensltvlty coefficient 

of dependent variable 

c 6 h 6 

c 6 h 5 

CO 

Pressure 


(a) T 0 = 1405K 

tp = 1 .0 Time = 220 ysec 


7 

CfjHt^ + Op ^ CfcHqO + 0 

-10.28 

-0.7590 

6.839 

0.1900 

10 

Ct^Hi; + Op C c;H^O + 0 

-.3129 

-.0680 

.4974 

.0220 

8 

C 6 H 5 0^ C 5 H 5 + CO 

-.2839 

-.0325 

.7193 

.0015 

6 

C 6 H 6 + OH C 6 H 5 + H 2 0 

-.7274 

.1342 

.3036 

.0053 

15 

c 6 h 5 ^c 4 h 3 * c 2 h 2 

-.0100 

-.0148 

.0449 

.0018 


(b) T 0 . 1600K 

tp = 2.0 Time = 140 ysec 


7 

C^Hg + 0 2 ^ C^H^O + 0 

-2.057 

-1.208 

2.663 

0.1927 

10 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 

-.1582 

-.1595 

.4286 


8 

C 6 H 5 0 ^ C 5 H 5 + CO 

-.0393 

■ 

.0938 


6 

C 6 H 6 + OH ^ C 6 H 5 + H 2 0 

-.3172 


.2488 

.0150 

15 

C 6 H 5 ^ C 4 H 3 + C 2 H 2 

-.2875 

-.6623 

1 .040 

.0930 


TABLE V. - SENSITVITY OF COMPUTED IGNITION DELAY TIME TO RATE CONSTANT 
VARIATION FOR TWO SHOCK IGNTIIONS 


Reaction 

number 

Reaction 

Tp, ys 

k std 

Tp, VS 

2 k std 



Percent 

change 

Tp, yS 

0.5 k std 

Percent 

change 


11 

0 

1 — 

to 

1405K, 

tp = 1.0 




7 

CftH^ + 0 2 ^ C^H^O + 0 

272 

188 

-30.9 

420 


10 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 


246 

-9.6 

313 

■ 

8 

C6H$0^ C 5 H 5 + CO 


268 

-1.5 

289 


6 

CfcHfc + OH ^ C 6 H 5 + H 2 0 


273 

.4 

280 

iBSSi 

15 

C 6 H 5 - C 4 H 3 + C 2 H 2 


272 

0 

276 

1.5 


(b) T 0 

= 1600K, cp = 2 

.0 



7 

C^Hc; + Op ^ C^H^O + 0 

151 

si 


203 

34.4 

10 

C 5 H 5 + 0 2 ^ C 5 H 5 O + 0 


IBB 


160 

6.0 

8 

C 6 H 5 0 ^ C 5 H 5 + CO 


Bn 

0 

154 

2.0 

6 

C 6 H 6 + OH ^ C 6 H 5 + H 2 0 


■B 

-1.3 

159 

5.3 

15 

c 6 h 5 -C 4 H 3 + C 2 H 2 


Mm 

-17.2 

172 

13.9 


TABLE VI. - BENZENE-OXYGEN-ARGON INITIAL CONDITIONS 



Equivalence 

ratio, 

4> 

Mole percent 

Initial temperature 
range, (behind 
reflected shock) , 

K 

c 6 h 6 

°2 

Ar 

1 

0.5 

na 

20.313 

78.333 

1209-1345 

2 

1.0 

Mm 

3.868 

95.616 

1345-1528 

3 

T .0 

119 

12.68 

85.63 

1283-1435 

4 

1 .0 

in 

12.582 

85.728 

1355-1408 

5 

2.0 

H9 

5.093 

93.555 
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TABLE VII. - COMPARISON OF COMPUTED AND EXPERIMENTAL IGNITION DELAYS 


Equivalence 

ratio 

<p 

Percent 

argon 

Initial 

temperature 

K 

Experimental 

Ignition 

delay, 

ysec 

Computed 

Ignition 

delay, 

ysec 

Percent 

differ- 

ence 

Standard 

deviation 

percent 


78.3 

1209 

878 

680 

-22.6 

52.5 




743 

600 

-19.2 




1254 

435 

498 

14.5 




1276 

330 

431 

30.6 




1291 

272 

388 

42.6 




1307 

185 

352 

90.3 




1314 

202 

334 

65.3 




1345 

159 

280 

76.1 


1 .0 

95.6 

1345 

755 

550 


22.0 

Dilute 


1374 

604 

450 





1402 

415 

370 





1412 

412 

345 

-16.3 




1428 

367 

313 

-14.7 




1482 

213 

216 

1.4 




1525 

122 

160 

31.1 




1528 

122 

159 

30.3 


1 .0 

85.6 

1283 

750 

600 

-20.0 

28.3 

Strong 


1290 

613 

575 

-6.2 




1294 

607 

556 

-8.4 




1328 

490 

450 

-8.2 




1355 

303 

378 

24.7 




1369 

287 

342 

19.2 




1379 

291 

320 

10.0 




1405 

198 

272 

37.4 




1408 

189 

266 

40.7 




1417 

178 

252 

41.6 




1435 

151 

225 

49.0 


2.0 

93.6 

1363 

1520 

870 

-42.8 

21 .6 



1415 

890 

607 

-31 .8 




1457 

599 

453 

-24.4 




1540 

274 

248 

-9.5 




1554 

243 

223 

-8.2 




1570 

211 

195 

-7.6 




1582 

157 

172 

9.6 




1600 

154 

151 

-1.9 



Standard deviation of all data points Is 33.2 percent. 


TABLE VIII. - COMPUTED AND EXPERIMENTAL 
ACTIVATION ENERGIES FOR IGNITION DELAY 


Equiv- 

alence 

ratio 

<p 

Activation energy cal/mole 

Percent 

differ- 

ence 

Experimental 

Computed 

0.5 

44160 

21230 

-52 

1.0 

Dilute 

41470 

27860 

-37 

1 .0 

Strong 

37250 

23560 

-37 

2.0 

42400 

32010 

-24 
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IGNITION DELAY, tlSEC PRESSURE, atm 



FIGURE 1. - COMPUTED PRESSURE AS FUNCTION OF TIME FOR BEN- 
ZENE-OXYGEN-ARGON SHOCK IGNITION. INITIAL TEMPERATURE, 
1363 K; EQUIVALENCE RATIO, <p, 2.0. 



1350 1280 1219 

TEMPERATURE, K 

FIGURE 2. - BENZENE-OXYGEN-ARGON IGNITION DELAY VERSUS RE- 
CIPROCAL OF TEMPERATURE FOR EQUIVALENCE RATIO 0.5. 



1515 1430 1350 

TEMPERATURE, K 

FIGURE 3. - BENZENE-OXYGEN-ARGON IGNITION DELAY VERSUS RE- 
CIPROCAL OF TEMPERATURE FOR EQUIVALENCE RATIO 1.0. 



1430 1350 1280 

TEMPERATURE, K 

FIGURE 4.- BENZENE-OXYGEN- ARGON IGNITION DELAY VERSUS RE- 
CIPROCAL OF TEMPERATURE FOR EQUIVALENCE RATIO 1.0. 


20 



IGNITION DELAY, msec 



1610 1540 1490 1440 1390 

TEMPERATURE, K 


FIGURE 5. - BENZENE-OXYGEN-ARGON IGNITION DELAY VERSUS 
RECIPROCAL OF TEMPERATURE FOR EQUIVALENCE RATIO 2.0. 



TIME, mSEC 

FIGURE 7. - TEMPERATURE VERSUS TIME FOR BENZENE-OXYGEN- 
NITROGEN REACTION. 
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MOLE FRACTION x 10' 



1515 1430 1350 

TEMPERATURE, K 


FIGURE 6. - BENZENE-OXYGEN-ARGON IGNITION DELAY VERSUS RE- 
CIPROCAL OF TEMPERATURE: EFFECT OF ARGON DILUTION FOR 
EQU I VALANCE RATIO 1.0. 



0 20 40 60 80 100 120 


TIME, mSEC 

FIGURE 8. - CYCLOPENTAD I ENE MOLE FRACTION VERSUS TIME FOR 
BENZENE-OXYGEN-NITROGEN REACTION. 



MOLE FRACTION X 10* MOLE FRACTION 



0 20 40 60 80 100 120 


TIME/ msEc 

FIGURE 9. - PHENOL MOLE FRACTION VERSUS TIME FOR BENZENE- 
OXYGEN-NITROGEN REACTION. 



TIME, msec 


FIGURE 11. - CARBON MONOXIDE MOLE FRACTION VERSUS TIME FOR 
BENZENE-OXYGEN-NITROGEN REACTION. 



0 20 40 60 80 100 


TIME. msEC 

FIGURE 10. - BENZENE MOLE FRACTION VERSUS TIME FOR BENZENE- 
OXYGEN-NITROGEN REACTION. 



TIME, msEC 


FIGURE 12. - CARBON DIOXIDE MOLE FRACTION VERSUS TIME FOR 
BENZENE-OXYGEN-NITROGEN REACTION. 
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